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ABSTRACT 

We follow a formalism presented by Kaiser [1] to calculate bulk flows in large 
scale surveys. We apply the formalism to a mock survey of Abell clusters a la Lauer 
& Postman [2] and find the bulk velocities in a universe with CDM, MDM and 
IRAS-QDOT power spectra. We calculate the velocity variance as a function of the 
f-D velocity dispersion of the clusters and the size of the survey. 

1. Introduction 

Recently Lauer and Postman (LP) [2] presented an analysis where they measured the reflex 
motion of the local group with respect to the 15000 km/s Abell Cluster frame. The velocity 
inferred is a 3<r result compared to the local group velocity in the CMB frame. They conclude 
that if the CMB dipole is doppler in origin one can infer the bulk velocity of the Abell Cluster 
frame in the cosmic frame to be 730 km/s. They made a peculiar velocity survey using the 
brightest cluster galaxy (BCG) of 120 Abel clusters as distance indicators with peculiar velocity 
dispersion of 16% of the radial redshift distance. Their survey is volume limited and complete 
to all clusters with an elliptical BCG in a sphere of 15000 km/s radius. They found that the 
L — a relation showed a dipole in the sky directed 80° from the direction of the CMB dipole 
and contend that this velocity of the Abell cluster frame relative to the cosmic rest frame (as 
defined by the CMBR) represents the true velocity on these scales. Given present theories it is 
quite surprising to find as large bulk velocities on such scales. Strauss [3] presents Monte-Carlo 
simulations of peculiar velocity data to find the probability that different theoretical models 
can explain the data, we approach the problem semi-analytically to test various cosmological 
models and survey results. 

We present an analysis based on work done by Kaiser [1] to infer what amplitude bulk 
flows one should measure given a specific velocity power spectrum. We show results for power 
spectra from the IRAS-QDOT survey [4], the BBKS CDM model [5] and from simulations of 
mixed dark matter (MDM) [6]. We show that the velocities are quadrature superposition of 
the noise that is inversely proportional to the number of data points and a convolution of the 
velocity power spectrum and the window function which is determined by the geometry of the 
survey. Clearly, the more power there is on large scales, the larger the expected velocities one 
would observe, therefore we chose as one of the power spectra we analyze the IRAS-QDOT 
plus the la error bars. 

In §2 we present the analysis and calculate the window function for an ensemble of observers 
surveying a similar sample to the LP survey and present the power spectra we used. In §3 we 
discuss the large N limit analytically. In §4 we present the results for a LP type survey, extend 
the analysis to similar surveys of different sizes and accuracy and compare with analytical 
results. We conclude in §5. 



2. Analysis 



We would like to calculate what one would expect to measure for the net streaming velocity 
of a sample of clusters given a universe with a velocity power spectrum P v (k). To this end, 
we assume an ensemble of observers each of whom measures the line-of-sight velocities of each 
cluster in a "mock" cluster sample constructed to mimic the features of the actual sample 
observed. 

Mock volume limited survey 





Fig. 1. — Aitoff-Hammer projection of a sample survey 



For comparison with the LP survey, we construct a volume limited survey by placing 
N = 120 clusters at random within a sphere of radius R = 150h _1 Mpc. (Below we will study 
how our results depend on the value of R.) As in the actual survey, we eliminate any cluster 
with galactic latitude < 15°. A typical "mock" survey constructed in this way is shown in Fig. 
1. 

The observed line-of-sight velocities for our sample clusters are assumed to have two sources 
of uncertainty; we assume that the peculiar velocities of the clusters are Gaussian distributed 
about the uniform streaming motion with 1— D dispersion <r* and that the observational errors 
are also Gaussian distributed and given by o n for the nth cluster. We shall take o n = f3r } so 
that the error in measurement is a fraction f3 / H of the distance r to the cluster, <r* to be a 
constant over the sample and H = lOOkm/s/Mpc. 

Given these assumptions, it can be shown [1] that the maximum likelihood solution for the 
uniform streaming motion Ui of a cluster sample is given by 



s 



where 



Here the cluster labeled by an index s has position f s and estimated line-of-sight peculiar 
velocity S s that is related to the true velocity by 



where e s is drawn from a Guassian with zero mean and variance (T 2 e = a 2 + a 2 . 
Writing the uniform bulk velocity in terms of Eqs. (1), (2) and (3) we get 
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Since the two terms are statistically independent, we can write the covariance matrix as 

Rij =< up. >= ^(J) + Eif . 
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Fig. 2. — normalized diagonal elements of the squared tensor window function for the survey 
and a Gaussian oc e~ k2R * 



The velocity term can be written as 



U\ v) (r ) = J d 3 r W %3 (r) v 3 (r + r) , 
where Wij is the tensor window function for the survey, given in fourier space by 



(6) 



W tJ Ck) = A~l e^. (7) 

s G s ~r ^* 



Then for 



wfm = w tm (k) w jn (k) k m k n 



the velocity part of the covariance matrix is a convolution of Eq. (8) and the velocity power 
spectrum 

i$> = / d 3 k Wf 3 (k) P v (k) , (9) 
where the velocity power spectrum is 

P v (k)=<\v(k)\ 2 >=^-P(k) . (10) 

We use the density power spectrum P(k) from the IRAS-QDOT survey [4] and its la result; 
the unbiased CDM transfer function (BBKS) [5] (Qh = 0.5, a 8 = 1); and the MDM simulations 
[6] normalized to the COBE quadrupole Q 2 = 17//K which is equivalent to a biasing parameter 
of b = 1.5. 

In Fig. 2, we plot the window function calculated for an LP type sample constructed as 
described above. The window function is fit well by a Gaussian with 75h _1 Mpc. In Fig. 

3 we plot same window function superposed with power spectra from models of interest. 
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Fig. 3. — Velocity power spectra, the window function for the survey and a Gaussian oc e * 
The noise term is 

Thus the final covariance matrix is 

R tJ = JtfkWlfaPvfa + Aj (12) 

The quantity 

A^^~; K^) = 4W\ (13) 
gives the variance of the distribution for the magnitudes of uniform streaming velocities for a 
given sample. Given A one can calculate the probability for observing a streaming velocity of 
a given magnitude for a particular survey type and power spectrum. 



3. Large N Limit 



In the limit ol large iV, we can replace many of the sums in section 2 with integrals. This 
gives us an alternative method of calculating the covariance matrix which acts as a check on 
the analysis described above. In what follows, we shall assume that we have N ^> 1 clusters 
distributed randomly in a sphere of radius R. For simplicity we shall concentrate on the two 
cases <t* = and cr s = 0. For the first case we shall take cr s = fir. 

To lowest order in 1/iV, we can write Eq. (2) as 

A; ^ + O(^) ■ (14) 

where we have used the fact that in the large N limit, A should be isotropic so that Aij = 
^ij A nn /3. For our two cases of interest 
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Next we will calculate Wij(k) using this method. Here we shall assume that the power 
spectrum P v (k) is a function only of the magnitude of k. With this assumption, we need 
only calculate the angle averaged Wij(k). In what follows we shall take all quantities to be 
averaged over angles and to be functions only of the magnitude of k. In this approximation 
Wij(k) = 6 tl Y. n W nn (k)/3. We then write Eq. (7) as 
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where we have taken k = kz for simplicity of calculation. 

In our angle averaged approximation, we make the replacement k m k n = 6 mn /3 in Eq. (8). 
Performing the integrals and using the expressions for A^ derived above we then find 



? (h] _ J % 3 (sm(kR) - kRcos(kR)Y /{kRf for a m = 
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where Si(x) = Jq dy sin(y)/y. 

Given a power spectrum, we can use the A^ and W^-(&) derived above to calculate 
and thus A in the large N limit. 

In realistic surveys, full sky coverage is never achieved. For the LP survey, for example, 
a 15° wedge about the galactic disk is removed, thus covering only 75% of the sky. For 
the noise term, A" 1 , we can account for this by multiplying our result by l/sqrt(0.75). For 
the other term, things are somewhat more complicated. We shall explore this question more 
thoroughly in an upcoming paper [7]. 



4. Results 



We have applied the formalism of section 2 to a set of "mock" surveys selected to mimic 
the features of the LP survey. For a given mock survey and power spectrum, the variance of 
the streaming velocity A can be calculated as described above. Due to statistical fluctuations, 
the value of A thus calculated will vary for different random realizations of a survey. From §3 
above we expect the variation in A over different realizations to be of order A/y/N ~ O.f A. 




200 300 400 500 600 200 300 400 500 600 
A (km/s) A (km/s) 



Fig. 4. — The distribution of A for the four power spectra we used for both (the true 
variance) and A (the left and right distribution respectively). The parameters for the Gaussian 
fit are shown. 

In Fig. 4 we show the result of 300 realizations of a LP type survey, with (3 = I6km/s/Mpc, 
cr* = 300km/s, radius R = I5000km/s and N = 120 clusters. The distribution of A over 
the different realizations is well fit by a Gaussian e -( A - A ™) 2 / 2<j2 anc [ j n good agreement with 
expectations. In Fig. 5 we show the variances as a function of the radius of the survey ball and 
the percent uncertainty (3 for the four velocity power spectra. As expected, the true variance 
(i.e. (3=0), indeed decreases as the size of the frame increases. However, for more realistic 
surveys, where the effect of the noise is hard to estimate accurately, the variance increases at 
large radii. 

In Fig. 6 we compare the results from the numerical analysis to the analytical results 
presented in §3 for the CDM[5] power spectrum with fi/j = 0.5 and <r 8 = I. As can be seen, 
the numerical and analytical analyses give similar results. The discrepancy for small R arises 
because in the analytical analysis we set <r* = for cr s ^ which has an affect in the small R 
regime. 

In the table below we give the results for the measured variance A and the noise-free 
variance A^. Next to them we present the likelihood of getting the LP velocity of 730km/s 
and the chance of getting this result, i.e. 100— CL (confidence level) of the result. The chances 
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Fig. 5. — The expected values of A for the four power spectra we used as a function of the 
maximum radius of the survey and the percent error in the the velocity dispersion. The solid 
curve is the true variance A^, the dotted, short dashed and long dashed lines are for {3 / H = 
10%, 20% and 30% respectively. 

of finding large amplitude velocities in a survey of the LP type are quite reasonable > 10% (i.e. 
these are 1.3 — 1.5a results) given the power spectra we used. However, the chance of these 
velocities to be the actual, true velocity for the frame are quite small < 2% (2.3 — 3.9a result). 
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Fig. 6. — The numerical and analytical results of A for the CDM power spectrum (fi/j = 
0.5 a 8 = 1) used as a function of the maximum radius of the survey and the percent error in 
the the velocity dispersion f3. The solid curves are the numerical reults, the dashed line are the 
analytical ones. 

5. Conclusions 

We have presented a formalism to calculate the velocity variance as a function of the 
geometry of the survey. We show that the variance is a quadrature superposition of the true 
velocity and the noise inherent in the observation. We apply the formalism to realizations of 
surveys similar to the Lauer & Postman analysis [2]. We show that comparing there results 
to the expected variance for this type of survey we find that there velocity is a 1 — 2a result 
depending on the power spectra we used, whereas comparing it to true expected velocity we 
get a 2.3 — 3.9cr result. It is indeed true that we do not know the power spectrum for this type 
of clusters and thus are unable to make a definitive statement about the expected velocities in 
this kind of frame. However, we believe that it is not unreasonable to assume that the power 
spectra we used are quite good estimates, and the true power is not significantly larger on 
scales > 75h _1 Mpc than the la IRAS-QDOT one. We therefore conclude that although it is 
not unlikely that LP would measure a velocity of the magnitude they quote without the need 
for extra power on large scales, the true velocity of the frame relative to the cosmic rest frame 
is likely to be significantly lower. Clearly, if the LP results are confirmed and represent the 
actual bulk flow of the Abell frame then the models presented here will be ruled out at high 
confidence level; further, any susccessful model will have to address the discrepency between 
the IRAS-QDOT and LP results. We will present a more rigorous analysis of the problem and 
include the clustering properties of the Abell catalogue in a subsequent paper [7]. 

Acknowlegement: We thank Nick Kaiser for valuable discussions and helpful conversations. 
HAF was supported in part by the National Science Foundation grant NSF-PHY-92-96020. 
RW was supported in part by NASA grant NAGW-2802. 



REFERENCES 



[1] Kaiser, N. 1988, MNRAS, 231 149-167 

[2] Lauer, T. Postman N. 1993, in this volume 

[3] Strauus, M. 1993, in this volume 

[4] Feldman, H. A., Kaiser, N. & Peacock J. A., 1993, UM AC 93-5, ApJ in press. 
[5] Bardeen, J. M., Bond, J. R., Kaiser, N. & Szalay, A. S. 1986, ApJ 304 15-61 
[6] Klypin, A., Holtzman, J., Primack, J. & Regos, E. 1992 SCIPP preprint 
[7] Feldman, H. A. & Watkins, R., in preparation 



This preprint was prepared with the AAS IATj^X macros v3.0. 



